Is the low-£ microwave background cosmic? 
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The large-angle (low-£) correlations of the Cosmic Microwave Background exhibit several sta- 
tistically significant anomalies compared to the standard inflationary cosmology. We show that 
the quadrupole plane and the three octopole planes are far more aligned than previously thought 
(99.9% C.L.). Three of these planes are orthogonal to the ecliptic at 99.1% C.L., and the normals 
to these planes are aligned at 99.6% C.L. with the direction of the cosmological dipole and with the 
equinoxes. The remaining octopole plane is orthogonal to the supergalactic plane at 99.6% C.L. 
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Much effort is currently being devoted to examining 
the cosmic microwave background (CMB) temperature 
anisotropics measured vifith the Wilkinson Microwave 
Anisotropy Probe (WMAP) [1-4] and other CMB ex- 
periments [5] . While the data is regarded as a dramatic 
confirmation of standard inflationary cosmology, anoma- 
lies exist. In particular the correlations at large angular 
separations, or low £, exhibit several peculiarities. 

Most prominent among the "low-£ anomalies" is the 
near vanishing of the two-point angular correlation func- 
tion C{9) at angular separations greater than about 
60 degrees. This was first measured using the Cosmic 
Background Explorer's Differential Microwave Radiome- 
ter (COBE-DMR) [6] and recently confirmed by obser- 
vations with WMAP [4]. This anomalous lack of large- 
angle correlation is connected to the low value of the 
quadrupole contribution, C2 , in a spherical harmonic ex- 
pansion of the CMB sky - 
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(with {2£ + l)Ce — J2m I'^^mP) ~ although the smallness 
of C2 does not fully account for the shape of C{6). The 
significance of the large-angle behaviour of C{6), espe- 
cially in light of the large cosmic variance in C2, is a 
matter of some controversy. The comparisons are, more- 
over, confused by the fact that where one author may 
calculate only the probability of the low value of C2 , oth- 
ers, such as the WMAP team, calculate the probability 
of C{d) being as low as it is over some range of angles. 
The issue of what prior probabilities and estimators [7-9] 
to use further complicates the statistical situation. 

While the overall absence of large-scale power has 
attracted the most attention, several other large-angle 
anomalies have been pointed out. De Oliveira-Costa et 
al. [11] have shown that the octopole is unusually pla- 
nar, meaning that the hot and cold spots of the oc- 
topolar anisotropics lie nearly in a plane. The same 
authors found that the axes r\2 and ns about which 
the "angular momentum" dispersion m^|a^m|^ of the 
quadrupole and octopole are maximized are unusually 



aligned, ]n2 • naj = 0.9838. Eriksen et al. [12] found 
that the deficit in large-scale power is due to a system- 
atic deficit in power between £ = 2 and 40 in the north 
ecliptic hemisphere compared to the south ecliptic hemi- 
sphere. Some of us [13] have shown that the £ = 4 to 8 
multipoles exhibit an odd, but very unlikely (~ 1% prob- 
ability), correlation with £ — 2 and I ~ i. These low-£ 
anomalies (and others [14]) have all been pointed out be- 
fore, but no simple connection has been made between 
them. Here we remedy that situation. 

By far the largest signal in the CMB anisotropy 
is the dipole, recently measured by WMAP [1] to be 
(3.346 ± 0.017)mK in the direction {I = 263°85 ± 0?1, 6 = 
48°25 ± 0°04) in galactic coordinates. This is caused by 
the motion of the Sun with respect to the rest frame de- 
fined by the CMB. As shown by Peebles and Wilkinson 
[15], the dipole induced by a velocity v is T{v/c)cos9, 
where is measured from the direction of motion. Given 
f = (2.725 ± 0.002)K [16], one infers that v ~ 370km/s. 

The solar motion also implies [15, 17, 18] the presence 
of a kinematically induced Doppler quadrupole (DQ). 
To second order in /3 = t;/c ~ 10~^, an observer mov- 
ing with respect to the CMB rest-frame sees the usual 
monopole term with a black-body spectrum, a dipolar 
term oc P cos with a dipole spectrum and a quadrupo- 
lar term oc /3^(3 cos^ — 1) with a quadrupole spectrum. 
Higher multipoles are induced only at higher order in 
f3 and so can be neglected. To first approximation the 
quadrupole spectrum differs very little from the dipole 
spectrum across the frequency range probed by WMAP. 
The DQ is itself a small contribution to the quadrupole. 
It has a total band-power of only 3.6/xK^ compared to 
(154 ± 70)fiK'^ from the cut-sky WMAP analysis [3], 
195. l^K^ extracted [11] from the WMAP Internal Linear 
Combination (ILC) full-sky map [19], or 201. 6^K^ from 
the Tegmark et al. full-sky map [10] (henceforward the 
Tegmark map). Therefore, it is a good approximation 
to treat the Doppler-quadrupole as having a dipole spec- 
trum plus a small spectral distortion which we shall ig- 
nore. We can then readily subtract the DQ from the ILC 
or Tegmark map. (The ILC and Tegmark maps differ in 
the amount of spatial filtering used to produce them.) 



2 




FIG. 1: The normal vectors for ^ = 2 (dark grey), and £ = 3 
(black), as well as ± the dipole direction (light grey), and 
the equinoxes (open circles) plotted in sinusoidal projection. 
The Galactic center is the coordinate origin. Galactic longi- 
tude is positive on the left. Grid lines are separated by 30°. 
The long-dashed line is the ecliptic; the short-dashed line is 
the supergalactic plane. The clustering of the i = 3 normal 
vectors is indicative of the "planarity of the octopole" [11]. 
(The clustering is clearer to the eye in the south hemisphere 
because of the projection.) Their closeness to the 1 = 2 nor- 
mal is indicative of the quadrupole-octopole correlation. Note 
how three normals are close to the dipole, the ecliptic and the 
equinoxes, while the fourth is on the supergalactic plane. 



Meanwhile, some of us showed [13] that the (.-\h mul- 
tipole, Ti, can instead be written uniquely in terms of 
a scalar A!^^^ which depends only on the total power in 
this multipole (i.e. on C^) and I unit vectors {v(^^*)|z = 
!,...,£}. These "multipole vectors" are entirely indepen- 
dent of C^, and instead encode all the information about 
the phase relationships of the a^m. Heuristically, 
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where v(^'*) is the i^^ multipole vector of the t"^ mul- 
tipole. (In fact the right hand side contains terms with 
"angular momentum" £—2, £—4, ... These are subtracted 
by taking the appropriate traceless symmetric combina- 
tion, as described in [13].) Note that the signs of all the 
vectors can be absorbed into the sign of AW. For conve- 
nience we take the vectors to point in the north galactic 
hemisphere. The multipole vectors for C. = 2 and 3 for 
the DQ-corrected Tegmark map are [in galactic {l,b)] 



^(2,1) ^ (11?26, 16?64), 

.(^(2,2) ^ (118'?87,25?13), 

v(3.i) = (22?63,9°18), 

■vi^'"^) = (86?94,39?30), 

v(3,3) ^ (-44?92,8?20). 
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(A similar analysis has been done for the ILC map. Re- 
sults are quoted where instructive.) 

As described in [13] there are several ways to compare 
the multipole vectors, however most striking is to com- 
pute for each t the ^{^—l)/2 independent cross-products. 



These are the oriented areas vvr(^'*'-') = ± (v(^'*) x v(^J)). 
The overall signs of the area vectors are again unphysical 
(we take them to point in the north galactic hemisphere), 
however their magnitudes are not. The area vectors for 
i = 2,2) for the Tegmark map (cf. Fig. 1) are 

,^(2,1,2) ^ 0.9900(-105?73,56?62), 

,^^(3.1,2) ^ 0.9017(-78?38,49?76), 

,^(3,2,3) ^ 0.9072(-141'?61,38?96), 

,^(3,3,1) ^ 0.9184(173?77,79°54). 
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The directions of w^^'^'^^ and of de Oliveira-Costa et aVs 
n2 are mathematically equivalent. The small difference 
is due to the removal of the DQ here in w^^'^'^^. 

Finally, the magnitudes of the dot products between 
^(2,1,2) ^jjjj yf(.3,i,j) ordered from largest to smallest are: 
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Using instead the normal vectors n(^'*'J) = 
w(^'''-'^/|w(^'*'-')|, the dot products are: 

Di= 0.9531, 1)2 = 0.8719, £'3 = 0.8377. (6) 

One can see from the large values of all three of the Dj 
that the quadrupole and octopole are aligned, and that 
the octopole is unusually planar. 

The Ai retain information about both the magnitudes 
and orientations of the w^^'*'-'^ . We have compared their 
values against 10'"' Monte Carlos (MCs) of Gaussian ran- 
dom statistically isotropic skies with pixel noise (as in 
[13]). The probability that the largest of these dot prod- 
ucts is at least Ai, the 2"^^ largest at least A2, and the 
3'''^ largest at least A3 is 0.021%. It is 0.11% without the 
DQ-correction, supporting that this is not just a statis- 
tical accident. (The ILC map yields 0.025% and 0.12%.) 

Ordering dot products does not induce a well-defined 
ordering relation for the MC maps [20]. A robust and 
more conservative statistic is the sum of the dot products, 
S = J2iAi. We find that only 0.128% of the MC maps 
have a larger S than the DQ-corrected Tegmark map. 
Thus the quadrupole-octopole correlation is excluded from 
being a chance occurrence in a gaussian random statis- 
tically isotropic sky at > 99.87% C.L. This result is sta- 
tistically independent of (though perhaps not physically 
unrelated to) the lack of power at large angular scales 
since all the information about the power is contained in 
the A'-'''\ and not in the multipole vectors. 

So far we have looked only at the correlation between 
the CMB quadrupole and octopole. Motivated by the 
results of Eriksen et al. [12], we next ask whether the 
quadrupole and octopole correlate with the ecliptic or 
the galaxy. We notice (see Fig. 1) that three of the four 
Yi(i-iJ) seem to lie near the ecliptic plane. Their dot prod- 
ucts with the north ecliptic pole are (for the Tegmark 




FIG. 2: Cosmic quadrupole plus octopole, TJ+s (from 
Tegmark map). Coordinates are as in Fig. 1, background 
color is 0/iK. The ecliptic (dashed line) avoids all extrema. 

map) - in order of size: 0.0271, 0.0450, 0.1786 and 0.5233. 
This means that three of the four planes defined by the 
quadrupole and octopole are nearly orthogonal to the 
ecliptic. The probability that a MC map has the same 
four dot products smaller than these is 0.104% (for the 
ILC map 0.193%). The sum of the four dot products, S, 
of MC maps exceeds that of the DQ-corrected Tegmark 
map only in 0.925%, so a chance alignment of the nor- 
mals with the ecliptic plane is excluded at > 99% C.L. 
A similar comparison to the galactic plane is unremark- 
able - the normal closest to the galctic pole, n^'^'^'^^ (the 
normal not near the ecliptic), could be even closer to the 
pole at 6.6% probablility; however, n'^^'^'^^ is only 0°07 off 
the supergalactic plane (see Fig. 1). Consequently, the 
likelihood of the dot product of one n^^'*'^^ with the su- 
pergalactic pole being as small as observed is just 0.449%. 

The three normals that are near the ecliptic also lie 
very near the axis of the dipole. The likelihood of this 
alignment with the dipole is 0.041% (0.098%) for the 
Tegmark (ILC) map and 0.405% with respect to the 
S statistics. Since the dipole axis lies so close to the 
equinoxes, this may also be viewed as an alignment with 
the equinoxes. The angular difference in ecliptic longi- 
tude of the normals from the equinox is < 1^, 13^^, 13'^, 
and 44^^. The probability of these four headless vectors 
having ecliptic longitudes this close to the equinoxes is 
0.039% (0.005%) for the Tegmark (ILC) map, 0.368% for 
the S statistics. This is statistically independent from 
their near orthogonality to the ecliptic poles. 

Two different statistical test (ordered dot products and 
sum of dot products) show evidence for low-^ anoma- 
lies in the Tegmark map (and consistent results for the 
ILC map). These anomalies could be due to a wrong 
treatment of foregrounds, although these full-sky maps 
do not refer to a foreground model. However, they em- 
ploy a linear combination of frequencies to produce a 
minimum variance map. They are susceptible to bias 
because chance alignments between CMB and Galactic 



3 




FIG. 3: Left to right: Cosmic (Doppler-subtracted) 
quadrupole, octopole, and T2+3 + T4 from the Tegmark map. 
The coordinate system and color scale are as in Fig. 2. The 
ecliptic (dashed line) is noticeably less correlated with these 
maps than with the quadrupole-octopole map of Fig. 2. 



signals tend to cancel in order to produce minimum vari- 
ance. Error estimates for the ILC and Tegmark maps 
have not been published, nevertheless we try to estimate 
the errors on the correlations induced by errors in the 
full-sky maps. As we add (ad hoc) Gaussian noise of 
10/iK to (which has a value of 9.3/iK in the Tegmark 
map), the evidence (from ordered dot products) for the 
quadrupole-octopole alignment is slightly reduced, from 
99.97% C.L. to 99.57% C.L. All other correlations re- 
main at > 99% C.L. Adding the same error to all aim 
with £ — 2 and 3 preserves all correlations at > 95% C.L. 

In Fig. 2, we plot a combined Doppler-subtracted 
quadrupole plus octopole map, T^^. The dashed line 
is the ecliptic. Strikingly, the ecliptic threads its way 
along the node line separating one of the hot spots from 
one of the cold spots, tracking the node over a third of 
the sky (and avoiding the extrema over the rest). The 
alignment is even better in the ILC map (not shown) 
than in the Tegmark map, which is not surprising given 
the spatial filter used by Tegmark et al. Replotting in 
other projections (not shown) does not noticeably alter 
the apparent correlation. A second look at Fig. 2 reveals 
a north-south ecliptic asymmetry - the three extrema 
in the north are visibly weaker than those in the south. 
One is cautioned that, given the observed planarity of the 
quadrupole-plus-octopole, one expects some hemispheric 
asymmetry because of the parity of the octopole. This 
may explain the asymmetry found in [12]. 

Interestingly, the correlation with the ecliptic is 
stronger in the combined quadrupole-octopole map, than 
in either separately, and the north-south ecliptic asym- 
metry is visible in neither alone. The correlation (but 
not the asymmetry) appears to be weakened if higher 
multipoles are added in, however a more thorough anal- 
ysis is merited. (For example, in the extrema seem to 
define both the ecliptic and the plane of the quadrupole- 
octopole, and the dipole direction lies near the center of 
a cold spot.) These trends are exhibited in Fig. 3. 

We have shown that the planes defined by the oc- 
topole are nearly aligned with the plane of the Doppler- 
subtracted quadrupole, that three of these planes are or- 
thogonal to the ecliptic plane, with normals aligned with 
the dipole (or the equinoxes), while the fourth plane is 
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perpendicular to the supergalactic plane. Each of these 
correlations is unlikely at > 99% C.L., and at least two 
of them are statistically independent. We have also seen 
that the ecliptic threads between a hot and a cold spot 
of the combined Doppler-subtracted-quadrupole and oc- 
topole map - following a node line across about 1/3 of 
the sky, and separating the three strong extrema from 
the three weak extrema of the map. 

We find it hard to believe that these correlations are 
just statistical fluctuations around standard inflationary 
cosmology 's prediction of statistically isotropic Gaussian 
random aemS- That the quadrupole-octopole correlation 
just happened to increase by ^ 5 when the quadrupole 
was Doppler-corrected seems particularly unlikely. The 
correlation of the normals with the ecliptic poles suggest 
an unknown source or sink of CMB radiation or an un- 
recognized systematic. If it is a physical source or sink 
in the inner-solar system it would cause an annual mod- 
ulation in the time-ordered data. An outer solar-system 
origin (beyond 100 A.U.) might avoid such a signal. It 
seems likely that a local source (sink) would also show up 
in polarization maps, especially there would be no reason 
for B-modes being significantly suppressed with respect 
to E-modes. Physical correlation of the CMB with the 
equinoxes is difficult to imagine, since the WMAP satel- 
lite has no knowledge of the inclination of the Earth's 
spin axis. Alternatively, this normals-ecliptic correlation 
could be a consequence of the closeness of the dipole to 
the ecliptic plane and the correlation of the normals to 
the dipole. But since the dipole is itself believed to be 
due to local motion, this would suggest non-cosmological 
contamination oi £ — 2 and 3, as does the observed cor- 
relation with the supergalactic plane. 

Although full-sky maps such as the ones we have used 
are not expected to be reliable at high 1 and do not agree 
with each other with respect to the galaxy, for the ex- 
isting maps (ILC map, Tegmark map and Eriksen et al. 
map [8]) the ecliptic plane does not show up in the dif- 
ference maps and their power spectra are consistent with 
each other at low-1 [8]. 

If indeed the £ = 2 and 3 CAIB fluctuations are incon- 
sistent with the predictions of standard cosmology, then 
one must reconsider all CMB results within the stan- 
dard paradigm which rely on low ^'s, including: the high 
temperature-polarization correlation Cf^ measured by 
WMAP [1] at very low £ (and hence the inferred redshift 
of reionization) ; the normalization of the primordial fluc- 
tuations (which relies on the extraction of the optical 
depth T from low £ ); and the running dng/dlogk of the 
spectral index of scalar perturbations (which, as noted in 
[21], depends on the absence of low-£ TT power). 
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